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Abstract. Differentially-expressed genes (DEGs) correlated 
to oral contraceptives (OCs) were identified by comparing 
the transcriptomes of cervical cancer patients who have taken 
OCs and those who have not. Their biological functions 
and relevance to clinical manifestations were investigated 
further in order to gain an understanding of the pathogenesis 
of cervical cancer and provide potential therapeutic targets. 
Level 3 RN A- sequencing (seq) data for cervical squamous-cell 
carcinoma and endocervical adenocarcinoma and the clinical 
information were downloaded from The Cancer Genome Atlas. 
The present study analyzed the RNA-seq data and information 
on OC use of 35 patients [OC users (n=18) and those who have 
never used OCs (n=17)]. Student's t-test was used in order to 
identify DEGs and the false discovery rate (FDR) was esti- 
mated by a Beta-Uniform Mixture model, which was adopted 
in multiple testing corrections. A functional enrichment 
analysis was performed with the Database for Annotation, 
Visualization and Integrated Discovery tool and BioCarta. A 
total of 80 DEGs were identified in OC users while FDR=0.3 
was set as the cut-off value. The metabolic process and human 
telomerase RNA gene transcription were significantly upregu- 
lated in DEGs. Furthermore, secreted LY6/PLAUR domain 
containing 1 was identified to be correlated to the patho- 
logical response, while the synapse defective 1 Rho GTPase 
homolog 2 was found to be significantly associated with the 
histological grade and overall survival time. In conclusion, 
present study shed light on the effect of OC use on the onco- 
genesis of the cervix and may indicate novel approaches for a 
targeted therapy of cervical cancer. 
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Introduction 

Cervical cancer is the second most common type of cancer in 
females worldwide (1). In China, there is an increasing trend 
in the incidence and mortality of cervical cancer. Numerous 
risk factors have been linked to this disease, including sexu- 
ally transmitted diseases, sexual debut time, menstruation 
and childbirth, smoking, immune deficiency and viral infec- 
tion (2). 

Oral contraceptives (OCs) are one of the widely used 
contraceptive methods for females. As the main components 
of OCs are hormones (estrogen and progesterone), OC use 
may alter the gene expression and thereby increase the 
incidence of cervical cancer (3). Analysis of data from a 
total of 28 studies including 12,531 females with cervical 
cancer indicated that the risk of cervical cancer decreases 
following termination of OC use (4). Another pooled study 
including 16,573 females with and 35,509 with no cervical 
cancer reported that the relative risk of cervical cancer was 
increased in current users of OCs and declined after their use 
was terminated (5). 

The mechanism underlying the increased risk of cervical 
cancer in OC users remains elusive. However, the most 
prevalent mechanism is a correlation between OC and 
human papilloma virus (HPV). HPV infection is thought to 
be the main but not the sole cause of the majority of cases of 
cervical cancer. Steroid contraception has been postulated 
to be one mechanism by which HPV exerts its tumorigenic 
effect on cervical tissue. It is believed that steroids bind to 
specific DNA sequences within transcriptional regulatory 
regions on the HPV DNA, either to increase or suppress the 
transcription of various genes (6). In females using combined 
OCs, the cervical mucus remains scanty, thick and highly 
viscous (7). Therefore, Guven et al (8) hypothesized that 
cervical mucus changes may contribute to the pathogenesis 
of cervical cancer, since scanty, thick and highly viscous 
cervical mucus may modulate and prolong the effect of 
carcinogenic agents. 

To assess changes in global gene expression resulting from 
OC use, the transcriptome of cervical cancer patients with and 
without OC use were compared and differentially-expressed 
genes (DEGs) correlated to OC use were identified. Further 
investigations were performed on their biological functions 
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and their relevance to clinical manifestations, which may be 
beneficial for the development of novel therapies. 

Materials and methods 

RNA-seq data and clinical information. Level 3 RNA-seq data, 
which had been obtained using a Illumina HiSeq (Illumina, 
San Diego, CA, USA), for cervical squamous-cell carcinoma and 
endocervical adenocarcinoma as well as the clinical informa- 
tion were downloaded from The Cancer Genome Atlas (TCGA; 
National Human Genome Research Institute, National Institute 
of Health, Bethesda, MD, USA; http://cancergenome.nih.gov/). 
The raw data provided reads per kilobase of exon per million 
mapped reads (RPKM, a measure of relative gene expression) 
for each gene. Data for 35 samples were obtained (17 samples 
from patients who never used OC, 16 samples from former users 
and 2 samples from current users). Given the small sample size, 
samples from former and present users were combined. Clinical 
data were also collected, including the time of diagnosis, time of 
death and follow-up time, from which the survival time could be 
calculated. The pathological response status, histological grade, 
age and race were also contained in the Biotab file. 

Screening of DEGs. R software, version 3.0.1 (9) was used 
for statistical analysis. A t-test was used to identify DEGs. 
In multiple testing corrections, the false discovery rates 
(FDRs) were estimated using a Beta-Uniform Mixture (BUM) 
model (10). 

Functional enrichment analysis. Gene Ontology (GO) enrich- 
ment analysis was performed with the Database for Annotation, 
Visualization and Integrated Discovery tool (DAVID; National 
Institute of Health; http://david.abcc.ncifcrf.gov/) (11) with 
hypergeometric distribution. The P-value was adjusted by the 
Benjamini & Hochberg method (12). BioCarta (San Diego, 
CA, USA; http://www.biocarta.com) was also used to identify 
enriched regulatory networks in DEGs. All statistical analyses 
were conducted using R software, version 3.0.1 (9) 

Gene expression and clinical manifestations. The relevance 
of gene expression to clinical manifestations was investigated. 
The t-test was applied to examine the correlation between 
genes and the pathological response. Analysis of variance 
(ANOVA) was used to analyze the gene expression among 
different histological grades. The Cox model was selected 
to explore the significance of the correlation between gene 
expression and the survival time. 

Results 

Differentially -expressed genes. The distribution of RPKM was 
similar for each sample following data pre-treatment (Fig. 1). 
DEGs were screened using the t-test. The distribution of the 
P-value and the fitted BUM curve are shown in Fig. 2. A total 
of 80 DEGs were disclosed, while FDR=0.3 was set as the 
threshold value. The fold changes of these genes are shown in 
Fig. 3. 

Cluster analysis results. A cluster analysis was performed using 
DEGs. The distance was defined as 1-Pearson correlation coeffi- 
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Figure 1. Box plot for the RNA-seq data. RPKM, reads per kilobase of exon 
per million mapped reads. 
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Figure 2. Distribution of P-values and fitted BUM curve. BUM, Beta-Uniform 
Mixture. 



cient and the Wards method was applied in hierarchical clustering. 
The result is shown in Fig. 4. A total of three samples from former 
OC users were grouped with the control. It is speculated that the 
disturbance of gene expression by OCs may diminish with time. 

Enriched biological functions. A GO enrichment analysis 
indicated that cellular macromolecule metabolic processes 
were over-represented in DEGs (Table I). BioCarta analysis 
indicated that DEGs were involved in the transcriptional regu- 
lation of the human telomerase RNA gene (hTERC) (Fig. 5; 
P=0.038, Fisher's exact test). 

Relevance of secreted LY6/PLAUR domain containing 1 
(SLURP 1) and synapse defective 1 Rho GTPase homolog 2 
(SYDE2) to clinical manifestations. SLURP 1 was the most 
upregulated gene, whereas SYDE2 was the second most down- 
regulated gene. Their relevance to clinical manifestations was 
explored. SLURP 1 was correlated to the pathological response 
(Fig. 6; P=0.014, t-test), while SYDE2 was associated with 
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Table I. Gene Ontology (GO) enrichment analysis result for the differentially-expressed genes. 
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T-value was adjusted with the Benjamini & Hochberg method. 
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Figure 3. Fold changes for each differentially-expressed gene. 



the histological grade (P=0.0028, ANOVA) and survival time 
(P=0.042, Log Rank test). An upregulation of SLURPl was 
observed in OC users and it was linked to the reduced frequency 
of the complete response (Fig. 6A). Low expression of SYDE2 
was observed in patients with histological grade 1 tumors 
(Fig. 6B) and it was linked with a longer survival time (Fig. 6C). 
Therefore, it was believed that SLURPl and SYDE2 were 
potential therapeutic targets for cervical cancer. Additionally, 
the inhibition of SYDE2 may improve the survival time. 

Discussion 

In the present study, the transcriptome of cervical cancer 
patients with and without OC use were compared in order to 



identify any relevant DEGs. A total of 80 genes were identified 
as DEGs, while FDR=0.3 was set as the cut-off value. This 
finding confirmed that OC use had a large impact on the gene 
expression in the cervix and thereby stimulated the pathogen- 
esis of cervical cancer. GO enrichment analysis revealed that 
the metabolic process was enriched in the DEGs. BioCarta 
analysis indicated that DEGs were involved in the hTERC 
transcriptional regulation. 

This finding was in accordance with previous studies, which 
have indicated that OCs have an impact on the lipid and carbo- 
hydrate metabolism (13). Hormones present in OCs significantly 
affect the plasma lipoprotein metabolism, which can raise the 
levels of plasma triglycerides, low-density lipoprotein, and 
high-density lipoprotein 3. They also affect the carbohydrate 
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Figure 4. Cluster analysis result. 



Downregulation of genes encoding for subunits of adaptor 
complex-3 has been reported in cervical cancer compared 
with normal controls (15). In the present study, AP3D1 was 
identified to be further downregulated in the OC group 
compared with non-OC group. AP3D1 is a subunit of the 
APS adaptor-like complex, which is associated with the golgi 
region and more peripheral structures. The AP-3 complex 
facilitates the budding of vesicles from the golgi membrane, 
and may be directly involved in trafficking to lysosomes (16). 
BMIl polycomb ring finger oncogene (BMIl) is a member of 
the polycomb group, which participates in axial patterning, 
hematopoiesis, cell cycle regulation and senescence. It may 
contribute to malignant cell transformation and its overexpres- 
sion has been reported in cervical cancer (17,18). However, 
according to the present study, BMIl was downregulated in 
OC users. Further studies on these genes may be helpful in 
elucidating the molecular mechanism. 

Telomerase is essential for the immortalization of the 
majority of human cancer cells. A study by Guilleret et al (19) 
indicated that hTERC expression is regulated during carcino- 
genesis. This gene is often amplified in certain types of solid 
human tumors, including the cervix (20). Guo etal (21) reported 
that hTERC amplification testing is a promising evaluation 
method for cervical cancer screening. This gene may be acti- 
vated by the transcription complex nuclear transcription factor 
(NF)-Y as well as the transcription factors specificity protein 



BIOCARTA 




Figure 5. hTERC transcriptional regulation was enriched in differen- 
tially-expressed genes according to the BioCarta analysis. hTERC, human 
telomerase RNA gene; pRB, retinoblastoma protein; NF-Y, nuclear transcrip- 
tion factor Y; SP, specificity protein. 



metabolism, primarily through the activity of progestin, causing 
conditions including insulin resistance, increases in plasma 
insulin levels and relative glucose intolerance (14). 

Further interesting DEGs were identified, including 
adaptor-related protein complex 3 delta 1 (AP3D1). 
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Figure 6. Relevance of SLURPl and SYDE2 to clinical manifestations. (A) Relevance of SLURPl to pathologic response (P=0.014, t-test); (B) relevance of 
SYDE2 to histological tumor grade (P=0.0028, analysis of variance) (Gl-4, grades 1-4); (C) relevance of SYDE2 to survival time (P=0.042, Log Rank test). 
SLURPl, secreted LY6/PLAUR domain containing 1; SYDE2, synapse defective 1 Rho GTPase homolog 2; KM, Kaplan-Meier. 



(Sp)l and retinoblastoma protein (pRB) and could be repressed 
by Sp3 (22). In the present study, Sp3 was downregulated in the 
OC group, which may contribute to the increased incidence 
of cervical cancer through pathways including hTERC. NF-Y 
alpha was also downregulated in OC users. Therefore, further 
studies are required in order to elucidate the effect of OC use 
on the hTERC transcriptional regulation. 

SLURPl and SYDE2 exhibited considerable fold changes 
in OC users and thereby, their relevance to clinical manifes- 
tations was investigated in the present study. The expression 
of SLURPl was identified to be negatively correlated with 
the pathological response, whereas the expression of SYDE2 
was negatively correlated with the survival time. SLURPl and 
S YDE2 may be potential therapeutic targets for cervical cancer. 

SLURPl is a member of the Ly6/uPAR family that 
lacks a GPLanchoring signal sequence. It is thought that 
this secreted protein contains antitumor activity. A study by 
Pettersson et al (23) indicated that SLURPl participates in 



the regulation of the gut immune functions and motility, as 
well as possibly having a role in colon carcinogenesis/cancer 
progression. Kalantari-Dehaghi et al (24) also identified that 
cancer-associated genes upregulation by nitrosamine 4-(meth 
ylnitrosamino)-l-(3-pyridyl)-l-butanone can be abolished 
by SLURPl (24). SYDE2, a GTPase activator, has also been 
reported to be a tumor suppressor (25). Methylated SYDEl was 
demonstrated to be significantly associated with the relapse-free 
survival of patients with breast cancer (26) However, there are 
currently no studies regarding the implications of SLURPl and 
SYDE2 in cervical cancer. Therefore, further studies should 
investigate their roles in this cervical cancer. 

In conclusion, the present study identified a range of 
DEGs associated with OC use in patients with cervical 
cancer. Further study of these genes may aid in elucidating 
the underlying mechanism of the generation of cervical cancer 
associated with OC use. Furthermore, SLURPl and SYDE2 
may be potential therapeutic targets for cervical cancer. 
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